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Abstract 

We study the large distance behavior of a steady distribution of two Brownian 
particles under external driving in a two-dimensional space. Employing a method of 
perturbative system reduction, we analyze a Fokker-Planck equation that describes 
the time evolution of the probability density for the two particles. The expression we 
obtain shows that there exists a long range correlation between the two particles, of 
type. 

1 Introduction 

The statistical properties of fluctuations at equilibrium are described by equilibrium statis- 
tical mechanics. This has been established through experimental measurements carried out 
to test the theoretical predictions of statistical mechanics. In contrast to the equilibrium 
case, there is no known general principle determining the statistical properties of fluctua- 
tions under nonequilibrium conditions. Indeed, it might be thought that it is quite difficult 
to obtain a universal theoretical framework on nonequilibrium fluctuations. 

There are many nonequilibrium steady states (NESSs) that settle down into an equilib- 
rium state if one condition, such as the strength of an external driving force or the chemical 
potential at a boundary, is controlled. In such NESSs, the statistical properties of fluc- 
tuations can be elucidated through an approach that seeks to determine how equilibrium 
fluctuations are modified under the infiuence of nonequilibrium conditions. 

In a pioneering work in this context, Kuramoto studied open chemical systems in 1974 
and pointed out that large scale fiuctuations should be considered separately from thermo- 
dynamic fiuctuations occurring locally in space pp. What is referred to in Ref. jTj as "long 
range coherence" is found more explicitly to appear in the form of long range correlations of 
fiuctuations for conserved quantities [2]. Here, there are two classes of systems that exhibit 
long range correlation. One class consists of systems driven by nonequilibrium boundary 
conditions. It includes laminar fiow systems jS], temperature gradient systems |3] and den- 
sity gradient systems |H] . In such systems, anomalous fiuctuations originate from the spatial 
inhomogeneity of averaged quantities. The power law decay exponents of the spatial correla- 
tion for conserved quantities are basically of l/r^'^"^) type, but different exponents can also 
be realized through composite effects |j6j. The other type of systems exhibiting long range 
correlation are locally driven systems. In such systems, the statistical properties of local 
fluctuations differ substantially from those of equilibrium systems. This modiflcation yields 



1 



long range correlation of the form in d (> 2) dimensional systems [7j. This behavior 
can be easily understood if we model the time development of a conserved quantity with a 
phenomenological linear Langevin equation |H1 IH] in which the anisotropy of both the current 
noise intensity and the transportation coefficient without detailed balance is assumed. 

In this paper, we inquire whether long range correlation of l/r'' type is peculiar to the 
fluctuations of macroscopic variables in driven systems. In general, a chain of correlated two- 
body interactions among many particles provides a contribution to the correlation function 
for the density fleld. Thus, when we consider the system of a microscopic level (at which 
particle motion is described), it is reasonable to conjecture that long range correlation ap- 
pears only in the macroscopic limit. However, if anisotropy with a local violation of detailed 
balance is the essence of long range correlation in driven systems, it may not be necessary 
to have a many body system in order to observe such correlation. As an extreme case, a 
system consisting of two particles under external driving may exhibit long range correlation. 

With this motivation, quite recently, in a calculation of the steady probability for the 
positions of two interacting random walkers in a d (> 2) dimensional lattice under external 
driving, it has been found that the large distance behavior of the probability, including the 
existence or non-existence of long range correlation, depends on the choice of the transition 
rules satisfying the condition of local detailed balance ^| . It is surprising that there is such 
a dependence, considering the fact that universal relations in the linear response regime 
do not depend on the precise nature of these rules. Given this situation, we are led to 
ask. Which rule is physically meaningful? However, it is difficult to answer this question 
by considering such transition rules themselves. For this reason, we investigate a physical 
model that corresponds directly to an experimental system. 

In the present paper, we study the large distance behavior of two Brownian particles with 
a local interaction under nonequilibrium conditions. Specifically, we consider two Brownian 
particles in a two dimensional space of temperature T. Their interaction is characterized by 
an interaction length ^. They are driven in one direction, which we choose as the x direction, 
by a constant external force /, under the influence of a periodic potential with period £. Let 
pi'\ri,r2) be the steady probability density for the positions of the two particles, (ri,r2). 
Calculating pi'\ri^ to leading order of in the asymptotic limit |ri — — oo under the 
assumption that ^ ^ £, we flnd that 

p(^)(ri,r2) -Ps(r-i)Ps(r-2) ^ ^^^^i^, (1) 

where c(ri, does not depend on |ri — r2|, but on the direction of the vector ri — r2, and 
Ps{t) is the one body steady probability density for a system without an interaction. Thus, 
it is found that a long range correlation of l/r"^ type appears in a system consisting of two 
locally interacting particles under external driving. 

In order to obtain this result, we apply a method of perturbative system reduction to the 
Fokker-Planck equation describing the time development of the probability density in the 
system described above. Here, what we refer to as "perturbative system reduction" consists 
of a perturbative calculation designed to obtain a simpler representation of a dynamical 
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system by restricting what we wish to describe. The first apphcation of perturbative system 
reduction to reaction diffusion systems was carried out by Kuramoto and Tsuzuki [T2]. They 
obtained a complex Ginzburg-Landau equation near a Hopf bifurcation and subsequently 
derived the simplest partial differential equation exhibiting spatially extended chaos [13j, 
which is now called the Kuramoto-Sivashinsky equation jXlj. These methods of derivation 
have matured since that time, and the universal structure underlying the calculations has 
been elucidated [HJ [T3j . 

This paper is organized as follows. In Section 2, we present the stochastic model we 
study and describe the basic features of the model. In Section 3, we formulate a perturbative 
system reduction by first reviewing the basic ideas introduced by Kuramoto. In Section 4, 
performing a perturbative expansion of the model, we obtain a reduced model describing the 
large scale behavior of the system in question. Using this result, in Section 5, we derive the 
asymptotic form of the large distance behavior of the steady probability density. Section 6 
is devoted to concluding remarks. 



2 Model 

We study the motion of two small particles (on the oder of micro-meters in radius) interacting 
with each other in a fluid of temperature T. The particles are confined to a two dimensional 
square of length L and are subject to a periodic potential U of period £ in a single direction, 
which we chose as the x direction. Typical systems with such properties can be realized 
experimentally fHl E|- Further, a flow with constant velocity can be used to apply a 
constant driving force / to the particles in the x direction. In this way, it is possible to 
experimentally realize NESSs for such a particle system. 

Let (T'i,r2) represent the positions of the particles. We assume that their motion is 
described by the Langevin equation 

where r = {x,y), V{ri — is an interaction potential with interaction length ^, (e.g. 
V{r) = const, for \r\ > ^.), and U{r) is a periodic potential satisfying 

U{r + ie,) = U{r). (3) 

Further, r}^{t) represents Gaussian white noise with zero mean and unit dispersion. Here, the 
Boltzmann constant is set to unity. For simplicity, we assume periodic boundary conditions 
in both directions and that the system size L is sufficiently larger than ^ and i. 

In this system, the probability density of the particle positions, p(ri,r2,t), obeys the 
Fokker-Planck equation 



dp 1 ^ d 

dt 7 ^ dvi 

1=1 



dU{ri) dV{ri - ra) dp 
-f^xP + P + P + T— 

OVi OVi OVi 



(4) 
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Although we wish to obtain the steady solution of it seems unfeasible to derive such a 
solution in exact form. For this reason, we formulate a perturbation method to extract the 
large scale behavior of the steady state solution under some assumptions. Before presenting 
the analysis, some preparation is needed. 

We first note that when V{ri — = 0, the steady state solution of (jlj can be derived 
easily in the form Ps{ri; f )ps{r2] /), where 

Psix,y;f) = ^I-{x), (5) 

with 

I_{x)= rrfx'e-'3^(^)+^^(^'+^)-'3^^'. (6) 
Jo 

Here, (3 = 1/T and Z is a normalization factor that is chosen so that we have 



dxps{x,yj) = i. (7) 







For later convenience, we define the following quantities: 



7 \ ox J "J ax 

Ur-J) = logPs(r;/). (9) 

Note that Js(/) does not depend on r, because Js(/) corresponds to the one particle steady 
state probability current for the case = 0. 

3 Formulation 

Next, we consider the effect of the interaction potential V by first writing 

p{ri, r2, t) = ps{ri; fi)ps{r2] f2)q{ri, ra, t), (10) 

where /i and /a are configuration dependent forces defined by 

N / dV{ri-r2) 
Mr,,r2) = f ^- . (11) 



Then, substituting (fTO|l into (jlj), we obtain the evolution equation of q{ri,r2,t) as 



^ = V 
dt ^ 

i=l 



Ps{ri;fi) dfi dxi 7% \dyi 



(12) 



where the operator Mj is defined as 



Ps{rufi)dxi 7 dxi dxi 'y \dx'f dyf ' ' 
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We study a simple situation where we fix the period ^ of the periodic potential U{r) 
and make the range ^ of the interaction longer and longer. We then introduce a small 
parameter e representing the extent of the scale separation: e = <^ 1. The existence of 
two typical length scales is account for explicitly by introducing the large scale coordinate 
Ri = {Xi,Yi) = evi and the periodic coordinate 9i = mod{xi,i), which expresses the i-th 
particle position in terms of the phase of the periodic potential. Using these coordinates, we 
rewrite U and V as 

U{r,) = U{d,) (14) 
Vir^-r2) = V{Ri-R2). (15) 

In a similar way, we further define Ps{di'i fi) and (psi^i', fi), where 

MR„R,)^f-e?^^^. (16) 

Now, we introduce a slowly varying field Q{Ri, R2,t) in such a way that 

q{ri, r2, t) = Q{R^, R2, t) + ep^'\e,, 62; [Q]) + ^V^'H^i, O2; [Q]) + ■ ■ ■ , (17) 

where 5'([Q]) represents the functional dependence of g on Q{Ri, R2,t). At this stage, Q is 
not determined. According to Ref. ^H] written by Kuramoto in 1989, Q can be regarded 
as the coordinate of a point on the slow manifold in the functional space {q}. With this 
interpretation, (fT7j) provides a representation of the slow manifold in terms of Q. Here, 
obviously, such a representation can be chosen arbitrarily. We therefore choose the following 
convenient form of the time evolution of Q is expressed by 

^ = eQ^'\[Q])+e'Q^'\[Q]) + ---. (18) 

If Q can be determined uniquely with this requirement, then ()18|) represents the system 
reduction we seek. Equations ()17|) and ()18|) constitute the basic assumptions of the per- 
turbative system reduction for ()12|) . with the replacements U by U, and so on. Thus, the 
problem we face is to determine whether p^'^\9i, O2] [Q]) and f2*^"'''([Q]) (n = 1, 2, • • • ) can be 
determined in an essentially unique way. In the next section, we see that indeed this can be 
done. 

4 Analysis 

We first substitute (jl7j) and (jl8|) into (jl2j) and extract all terms proportional to e. We then 
obtain 

2 

n^'\m = Y}Mf^p^'^ + M'^'^Ql (19) 

i=l 
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where the operators M^^^ and M^^^ are given by 



(0) 



Jsifi) d 



2 ' 



(1) 



■Ji) d T 



(20) 
(21) 



Because Ivlf^'^ -1 = 0, p^^^ can be obtained only when the solvabihty condition is satisfied. 
Thus, we need to find an exphcit form of the solvabihty condition. 
Let us define an operator Lj as 



Then, using the relation 

for an arbitrary square integrable periodic function (f){Oi), we obtain 

d9,pM){M^''Um) = 0. 



From this, the solvability condition for (fT^ turns out to be 



d9id92Ps{9i)p,{92 



which yields 



i=l 



0, 



i=l 



(22) 



(23) 



(24) 



(25) 



(26) 



Under this condition, we can derive p^^^ of (fT^ in the form 
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P 



(1) 



i=l 



dQ 
dX, 



X 



(27) 



where a{9i; fi) can be calculated explicitly as in Ref. (TH], and x^^K^i^ R2) is an arbitrary 
function of (i?i, -R2)- (The choice of this function does not infiuence the result for 
Next, we sum up all terms proportional to e"^. This yields 



+ 



^d(/.) d'v 
pMJ,)dXf 
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Q + 



1 d dV 



Q 



(28) 



where 



(2) 



+ 



7 V^^,' dYN ' 

and we have introduced the differential mobihty yUd(/) = dJs/df, which is found to be 



(29) 



X 



2 ' 



(30) 



with 



Ux] 



^^i^pU{x)-l3U{x-x')-Pfx' _ 



(31) 



The solvabihty condition for p^^-* in ()28p. which is obtained by multiplying both sides of 
by / <i#id«2P,(«i)p.(«2), yields 



Atd(/, 



9X 



Q + D{f) 



d^Q Td^Q 1 d dV 



+ 



where D is obtained as 



D{f) 



T\J^dx{l4x)fI, 



[X 
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(32) 



(33) 



This quantity clearly represents the diffusion constant in the x direction [TF. The derivation 
here is parallel to that given in Ref. Under the solvability condition (j32j) . we can obtain 

Finally, carrying out a similar calculation, we can determine p'^")(6'i, 62; [Q]) and ^^"•'([Q]) 
from the terms proportional to in (fT^ . with ()17|) and (fTBj) . This iterative procedure 
constitutes the perturbative system reduction. 



5 Long range correlation 

Recall that our goal is to obtain the steady state solution of the Fokker-Planck equation (jH) 
under the assumption e <^ 1. Let us express this solution as 



P^^\ri,r2) = Psiri)ps{r2)qsiri,r2). 



(34) 



Then, as far as we focus on the large distance behavior of gs(ri,r2), from (fTUj) . (fT7|) . (fTHj) . 
(j2ni) and we can assume that gs(ri,r2) satisfies the equation 



E 

j=i 



(9xf 7 %2 ^Qy. \gy. 



0. (35) 
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Now, we define ip{ri, r2), which represents the non-equihbrium contribution of Qs, through 
the relation 

Furthermore, assuming that V is sufficiently small, we linearize ()35p with respect to and 
V. This yields 



E 



1=1 



+ (/}(/) -/id(/)r) 



dxi 



dV dip 
-f3— + 



dxi dxi J 7 \dx^ dy, 



+ 



0. (37) 



Then, using the Fourier expansion, we can solve for i\} in (jHTj) as 



where 



V'(ri,r2) = p-/idT)/3 



V{k) 



-ifc(ri— r2) 



Dkl + Tklh 



V{k), 



rfVe*'=^V(r). 



From this result, it is straightforward to derive the asymptotic form 



1 D - /idT^i/2^^^^_ (xi - X2)yD - {y, - y^Yi/T 



27r (DT)3/2 



[[x,-X2f/D + {yr-y2fl/Ty 



(38) 



(39) 



(40) 



in the region ^ -C |ri — <^ As shown in Ref. jl9j, the Einstein relation D = n^T is 
violated in general NESSs. Thus, we conclude that there is long range correlation of 
type. 



6 Concluding remarks 

We have demonstrated that there exists long-range spatial correlation between two interact- 
ing Brownian particles under external driving. We have found that this long range correlation 
is proportional to D — fi^T, which represents the degree of the breakdown of detailed balance. 
It is quite reasonable to expect that the long range correlation found for this two particle 
system exists also in many particle systems, with a quantitative correction arising from many 
body effects. It is a future project to study many body effects by extending approach in the 
present paper. 

The existence of long range correlation of l/r'^ type makes it difficult to construct a 
universal framework for a statistical theory. Let us explain the reason for this by considering 

particles in a two dimensional box of length L. We write the A^-body steady state 
distribution for this system as 

Pi{r^}) = e-^*"'-^". (41) 

One may naively interpret $({rj}) as "effective energy " of the particles under nonequilib- 
rium conditions, because $({rj}) corresponds to the Hamiltonian of the particle configu- 
ration {rj} at equilibrium. With this interpretation, there appears an effective long range 
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interaction potential of type. Then, for such systems, simple considerations yield the 
estimation 

{^) r^L"^ log L (42) 

in the limit that L oo, with fixed L^/N. This implies that extensivity, which is the most 
essential property of thermodynamic systems, does not hold. Such statistical systems are 
pathological for the following reasons. First, a statistical distribution located in the central 
region of the system depends sensitively on the nature of the boundary conditions [20] • That 
is, it is difficult to define a bulk region for the system. Second, if $ can be measured as 
"energy" using some experimental method, ()42j) implies that a significant amount of energy 
can be extracted by merely splitting one system or combining two systems. Because such 
a situation seems to be unphysical, it is reasonable to conjecture that $ does not represent 
an "energy" . These conclusions cast doubt on the possibility of realizing a unified statistical 
framework. 

Despite these seemingly intractable properties, we wish to seek a universal statistical 
framework for NESS by separating the problem in the following way. First, we propose 
to check the possibility that large scale fiuctuations can be distinguished from small scale 
fluctuations in some way. If this can be done, we hope to determine whether small scale 
fluctuations, which may deviate substantially from those of an equilibrium system, can be 
characterized in terms of an energetic quantity. Then, flnally, we hope to study large scale 
fluctuations on the basis of the characterization of smalls scale fluctuations. Recently, we 
have made some progress in the characterization of small scale fluctuations through an 
extension of thermodynamic functions [211 122] • We now propose to attempt unifying large 
scale anomalous fluctuations with our thermodynamic framework, going beyond the result 
of the present study. 

In closing this paper, we would like to return to Ref. ^ written by Kuramoto in 1974. 
The observation made there that large scale fluctuations should be considered separately 
led Kuramoto to focus on dynamical behavior of macroscopic variables. In particular, when 
solutions of deterministic equations for macroscopic variables describe a rich variety of phe- 
nomena including oscillations and chaos, the understanding of such phenomena from a dy- 
namical system point of view may be most important. With this realization, Kuramoto 
naturally was led to study dynamical systems. This is regarded as the genesis of nonlinear 
dynamics as a method for studying nonequilibrium statistical phenomena. The most im- 
portant message here seems to be that to formulate questions that do not conform to the 
contemporary mainstream can lead to new flelds of research. Today, the study of nonlinear 
dynamics has been fully developed. Following Kuramoto, we should consider to seek the 
formulation of precise and deep questions regarding nonlinear and nonequilibrium systems 
that do not conform to current trends. 
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